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I. INTRODUCTION 



In the recent years, the intersecting D-brane configurations have played an important role 
in several areas. The most prominent one is the construction of four-dimensional solutions 
[3-7] of Type II string theory, compactified on orientifolds . In particular the appearance 
of the chiral matter [1, 2] at the brane intersection provides a promising starting point to 
construct models with potential particle physics implications [1, 2]. 

The model-building with intersecting branes was developed in a series of papers. In par- 
ticular, non-supersymmetric [3-7] (and subsequently explored in [8-11]) and more recently, 
super symmetric [12-18] constructions with quasi-realistic features of the Standard-like and 
grand-unified models have been given. One has a tremendous freedom in the construc- 
tions of non-supersymmetric models, since the Ramond-Ramond tadpole cancellation con- 
ditions can be satisfied for many brane configurations leading to the Standard-Model gauge 
group and three families of quarks and leptons. However, the fact that the theory is non- 
supersymmetric introduces the Neveu-Schwarz-Neveu-Schwarz uncancelled tadpoles as well 
as the radiative corrections of the string scale. (For the constructions with intersecting D6- 
branes the string scale is necessarily of the order of the Plank scale. However examples [19] 
with intersecting D5-branes have been given, where the string scale can be as low as the 
TeV scale.) 

On the other hand supersymmetric intersecting D-brane constructions are extremely con- 
straining. Nevertheless such supersymmetric constructions with intersecting D6-branes, 
which have the Standard-like [12, 13, 15, 17] and grand unified model spectra [13, 18], 
have recently been constructed. In particular, these models have an additional quasi-hidden 
gauge sector that is typically confining which may have interesting implications for the su- 
persymmetry breaking [22]. Note however, that these models typically suffer from additional 
exotics [20]. Both supersymmetric and non-supersymmetric constructions have adjoint mat- 
ter associated with each brane configuration, since the toroidal cycles wrapped by the branes 
are not rigid. Interestingly, the embedding of supersymmetric four-dimensional models with 
intersecting D6-branes has a lift [13, 14] into M-theory that corresponds to the compactifi- 
cation of M-theory on a singular G 2 holonomy manifold [12, 13, 25-27]. 

While phenomenology of both non-supersymmetric [6, 11, 23] and supersymmetric [20- 
22] models has been addressed, the actual string calculations of the couplings in this theory 
have been limited. While the tree-level gauge couplings are relatively easy to determine 
and their features have been studied, see e.g. [20] and references therein. A calculation of 
gauge coupling threshold corrections [24] is also of interest, since it could be compared to 
the strong coupling limit of M-theory compactified on the corresponding G 2 holonomy space 
[28]. 



2 



An important set of tree level calculations involves the open-sector states that appear 
at the brane intersections. These states include the chiral matter. In the supersymmetric 
constructions the appearance of the full massless chiral supermultiplet is ensured there. The 
couplings of most interest are the three linear superpotential couplings, such as the coupling 
of quarks and leptons to the Higgs fields. On the other hand the four-point couplings are 
also of interest, since they indicate the appearance of potentially other higher order terms 
in the effective Lagrangian. 

The calculations of couplings of states (tachyons) appearing at the non-supersymmetric 
intersections of branes also has interesting implications in the study of tachyon potential 
and the phenomenon of tachyon condensation. In particular for specific T-dual models 
of [p-(p+2)] bound state configurations the four-point calculations have been addressed in 
[31-33]. 

The purpose of this paper is to perform explicit string calculations of the four-point and 
three-point correlation functions associated with the states appearing at the intersections of 
branes that wrap cycles of the internal tori. The non-trivial part of the calculation involves 
the evaluation of the correlation functions of four (three) bosonic twist fields, which signify 
the fact that the states at the intersection arise from the sector with twisted boundary condi- 
tions on the bosonic (and fermionic) string degrees of freedom. (For supersymmetric cycles 
the physical massless states at the intersection correspond to the chiral supermultiplets.) We 
employ the techniques of conformal field theory, which are related to the study of bosonic 
twist fields of the closed string theory on orbifolds [29]. Similar techniques were employed 
in the study of Type II string theory for bound-states of p-(p+2) brane sectors [31-33]. 

Specifically, we focus on intersecting D-branes wrapping factorizable N-cycles of T 2N = 
T 2 x T 2 ■ ■ ■. Thus, in each T 2 the D-branes wrap one-cycles, and the problem reduces 
to a calculation of correlation functions of bosonic twist fields associated with the twisted 
sectors at intersections of D-branes on a general T 2 . Thus the final answer is a product of 
contributions from each correlation function on each T 2 . 

We provide a general result for: 

(a u (x 1 )a- u (x 2 )a l/ (x 3 )a- u (x 4: )) , (1) 

which corresponds to the bosonic twist field correlation function of states appearing at the 
intersection of two pairwise parallel branes with intersection angle 7rz/ (See Figure 1). 

In the case of the twist fields appearing at the same intersection, our result is interpreted 
in terms of the volume of the torus, the lengths of the one-cycles L\ and L 2 that each set 
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of branes wrap and the intersection numbers I. 1 We also address the case when the twist 
fields are associated with different intersections of the two branes. In particular we address 
in detail the summation over the instanton sectors for such general cases. 

The next calculation that we set out to do is that of the four-point correlation function: 



(a l/ (x 1 )a- l/ (x 2 )o-- X (x3)o-x(x4)), (2) 

which corresponds to the bosonic twist field correlation function of states appearing at the 
intersection of two branes intersecting at respective angles and n\ with the third set of 
parallel branes (See Figure 2). This correlation function is specifically suited for taking the 
limit of X2 — > £3 which factorizes to a three point function associated with the intersection 
of three branes. This latter result is particularly interesting since it provides a key element 
in the calculation of the Yukawa coupling. 

In this set of calculations we determine both the classical part and the quantum part 
of the amplitude and thus obtain the exact answer for the calculation. In particular, the 
calculation of the quantum part depends only on the angles (and is thus insensitive to the 
scales of the internal space). On the other hand the classical part carries information on the 
actual separation among the branes and the overall volume of T 2 as well. 

In particular the full expression (both classical and quantum part) for the Yukawa cou- 
plings for branes wrapping factorizable cycles of T 6 is written as 2 



Y = V2g 2ir T — — ; J ' ] > - - > exp 3 



where Aj(m) is the area of the triangle formed by the three intersecting branes on the 
j-th torus and go = e*/ 2 , with $ corresponding to the Type IIA dilaton. The coupling 
is between two fermion fields and a scalar field, i.e. the massles states appearing at the 
respective intersections, whose kinetic energies are taken to be canonically normalized. 

While we were in a process of completing this work the paper [34] appeared where a 
comprehensive analysis of the classical part of the string contribution (disk instantons) 
to the Yukawa coupling in models with intersecting branes on Calabi-Yau manifolds was 
given, and extensive explicit calculations of the classical string contributions for models of 
intersecting branes on toroidal orientifolds were presented. 



1 A special case, when the branes wrap the canonical cycles of a torus with the complex structure v, a 
T-dual interpretation of this correlation function is that of the bosonic twist fields for D0-D2 brane with 
the magnetic flux B = cot(7w). 

2 See Note Added at the end of the paper. 
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(a) 



(b) 



FIG. 1: Target space: the intersection of two parallel branes separated by respective distances d\ 
and c?2 and intersecting at angles irv (Figure la). World-sheet: a disk diagram of the four twist 
fields located at £1,2,3,4 (Figure lb). The calculation involves a map from the world-sheet to target 
space. 

Our work has certain overlap with that of [34]. In particular, our work focuses only 
on models with branes wrapping factorizable N-cycles of T 2N = T 2 x T 2 ■ ■ ■. We evalu- 
ate the classical action contribution by explicitly solving for the classical solutions of the 
bosonic string with the boundary conditions governed by the locations of the D-branes. For 
the special case of the three point function we therefore also derive the result of [34] that 
the classical string contribution to the three-point coupling involves a summation over the 
exp(— A/2na') : where ^-corresponds to the area of the triangles associated with the inter- 
sections of the branes in each T 2 . On the other hand, we have also determined the quantum 
part of the correlation functions, thus obtaining the full expression for the couplings. 

The paper is organized as follows. In Section II we determine the correlation function 
(1). In Section III we calculate (2) and factorize it on a three point function to determine 
the corresponding Yukawa coupling. Conclusions, that include comments on generalizations 
of these calculations as well as physical implications are given in Section IV. 
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FIG. 2: Target space: the intersection of two branes intersecting respectively with the two parallel 
branes at angles iris and tt\, respectively (Figure 2a). World-sheet: a disk diagram of the four 
twist fields located at £1,2,3,4 (Figure 2b). The calculation involves a map from the world-sheet to 
target space, allowing for a factorization on three-point function. 

II. FOUR-POINT FUNCTION WITH ONE INDEPENDENT ANGLE 

In order to evaluate the path integral for the partition function of open strings stretching 
between D-branes intersecting at an angle tiu we split the embedding fields X 1 = X l cl + X 1 
into a classical solution to the equation of motion, subject to the appropriate boundary 
conditions, and a quantum fluctuation. The mode expansion for the quantum fluctuation 
is not integer moded due to the boundary conditions. The vacuum of the X 1 CFT is then 
created by primary fields a v acting on the SL(2, R) invariant vacuum. The partition function 
naturally factorizes into a classical contribution due to worldsheet instanton sectors and a 
quantum amplitude due to quantum fluctuations. In contrast to the instanton contribution, 
the quantum amplitude contains no topological information about the worldsheet embedding 
in target space, but it is still essential for the complete determination of Yukawa couplings 
in a general model with intersecting branes. 

A. Evaluation of Quantum Amplitude 

We shall employ the stress tensor method [29, 31, 32] to evaluate the quantum amplitude 
of four twist operators. For oriented theories the twist operators live on the boundary of the 
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disk and change the boundary conditions as we move along the boundary. The boundary 
conditions are specified by the D-brane configuration in target space (See Fig. 1). We 
concentrate on a single T 2 and Dl-branes wrapping one-cycles. The amplitude for branes 
wrapping factorizable three-cycles on T 6 is then the product of the amplitudes for the three 
T 2 factors. 

In terms of the complexified coordinates X = X 1 + iX 2 , X = X 1 — iX 2 on T 2 the 
boundary conditions read 

dX + dX = 0, dX + dX = 0, on (-oo, a^) U (x 2 , x 3 ) U (x 4 , +oo) (4) 
e inu dX + e- i7TU dX = 0, e~ iwu dX + e iwu dX = 0, on (x 1 ,x 2 ) U (x 3 , x 4 ). 

These conditions define the OPEs of the embedding fields with the twist operators, namely 

dX(z)a u (x) ~ ( z - x y~ 1 T u (x) + ... (5) 
dX(z)a u (x) ~{z- x)- v t' v {x) + ... 
dX(z)a„(x) ~ -(z - x)- v t' v {x) + ... 
dX{z)a v {x) ~ -(z - xy- l r v {x) + ... 

and similarly for cr_„(x). To evaluate the correlation function of four twist fields 
(a u (x 1 )<j- u (x2)cr l ,(x 3 )(j- v (x4 : )) we consider the correlators 

(-^jdX(z)dX(w)a v (x 1 )a- v (x 2 )a v (x 3 )a- v (x 4 )) 

g{Z,W) = ; r r r rr (6) 

{a v (x 1 )a- v (x 2 )cr u (x 3 )a- u (x4)) 

_ (-^dX(z)dX(w)a u (x 1 )a^ l/ (x 2 )a l/ (x 3 )a^ u (x 4 )} 
(a l/ {xi)a_ u (x 2 )a u (x 3 )a_ u {x 4 )) 

(-^dX(z)dX(w)a l/ (x 1 )a_ l/ (x 2 )a l/ (x 3 )a_ l/ (x 4 )) 

g(Z,W) = ; -. r -. r -. r -. rr (8) 

{a v [x l )a^ v {x 2 )a v {x 3 )a_ v {x 4 )) 

and 

- ^ _ {--^dX{z)dX{w)a u (x l )a_ u (x 2 )a u (x 3 )a_ u (x 4 )) 
{<j v (xi ) G- v (x 2 ) o v (x 3 ) G- v (x 4 ) ) 
The OPEs (5) together with the conditions 

g(z, w) ~ (z — w)~ 2 , h(z, w) ~ regular (10) 

as z — > w uniquely determine 

g(z,w) = ui- v (z)u v (w) 



> x x (f ~ x i)( z ~ x ^)( w - x 2 ){w - x 4 ) | 
(z — w) 2 



(z - x 2 )(z - x 4 ){w - x^jw - x 3 ) | ' 

(z — W) 2 
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(a) (b) 

FIG. 3: World sheet contours. The contours C\ and C 2 (Figure 3a) are the two topologically 
inequivalent contours leading to two independent conditions. The contours in Figure 3b define the 
global monodromy conditions used in Section III. 

and 

h(z, w) — —u l/ (z)u l/ (w)B({x i }) (12) 

where 

u u (z) = {z- Xl )-»(z - x 2 y-\z - x 3 )-"(z - x 4 y-\ (13) 

Here A({ Xi }) and B({ Xi }) are functions of the twist field positions to be determined. The 
boundary conditions and holomorphicity imply 

h(z,w) = -g(z,w), g(z,w) = -h(z,w). (14) 

In order to determine the functions A and B we impose appropriate monodromy con- 
ditions which will ensure that the quantum fluctuations X qu are local. This is guaranteed 
if 

I dX = [ dX = (15) 

where Ci is any non-trivial worldsheet contour. In the case at hand there are two topolog- 
ically inequivalent contours C\ joining the intervals (£2,^3) an d (—oo,Xi) and C 2 joining 
the intervals (x 3 ,x 4 ) and (xi,x 2 ) (See Fig. 3a). One can save some effort, however, by 
noticing that the contours can be analytically continued to the worldsheet boundary along 
which only one particular linear combination of target space fields satisfies Neumann bound- 
ary conditions and can therefore have non-trivial displacement (along the boundary of the 
world sheet d = drd T and hence the fields satisfying Dirichlet boundary conditions give no 
contribution). The non-trivial conditions are 

/ d(e iwu X - e~™ v X) = (16) 

JCi 
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/ d(X-X) = 0. (17) 
Jc 2 

When inserted into the four twist field correlation functions these lead to monodromy con- 
ditions for the correlators g, g, h and h. For example, condition (17) implies 

/ dz[g(z,w) — g(z,w)] + / dz[h(z : w) — h(z, w)\ — 0. (18) 
Jc 2 Jc 2 

Now, using relations (14) we can trade h and h for g and g. Moreover, z can be traded for z 
by integrating along the mirror image of contour C2 about the real axis, call it contour C<i- 
Taking into account the phases of g and g on each of the contours one sees 

/ dzg(z,w)= / dzg(z,w), (19) 
/ dzg(z,w) = / dzg(z,w). (20) 

[ dz[g(z,w)-g(z,w)}=0. (21) 

J Co 



and hence 



Similarly we derive 

/ dz[e i7TU g(z, w) - e~ i7n/ g(z, w)} = 0. (22) 
Jc\ 

Invoking SL(2, M.) invariance we fix x\ = 0, X2 = x, x 3 = 1, x±^> oo. Dividing then by 
uj v {w) and letting w^oowe get 

g(z, w) -> [u(z -x) + A(a0]a>i_„(2) (23) 

g(z,w) ^B{x)Cj v {z) (24) 

where A(x) and _B(x) have been redefined appropriately to absorb the singularity arising 
from X4 — > oo and 

= (z - xO-^z - a*)"" 1 ^ - x 3 )- v . (25) 
Conditions (22) and (21) then give 

v \ dz(z — x)Cui- v (z) + A(x) / dzuj\- v {z) — B{x) I dzuo v {z) = (26) 

J Co. J Co J Co 



and 



v I dz{z - xfa-viz) + A(x) dzG)x- v (z) -e~ 2wiu B(x) dzu„(z) = 0. (27) 
J a Jc^ Jc^ 



'Ci JC*i JCi 

Evaluating the contour integrals we find 



—v [ dz(z — x)Q\- v {z) — x(l — x)-j- [ dzQ\- v (z), (28) 
Jc, dx J c . 
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/ dzu^z) = F(x) (29) 
Jc\ 

and 

dzui- v (z) = e i7TU F(l - x) (30) 



where 



F(x) = B{v,\-v)F(y,\-v',\;x)= [ dyy'^l - y y-\l - xy)~ v . (31) 

Jo 

B(u, 1 — v) is the Euler Beta function and F(a, b; c; x) is the Hypergeometric function. 
Solving for A(x) we finally obtain 

A(x) = l -x{\ - x)d x \og[F{x)F{\ - x)]. (32) 

The quantum contribution Z qu (x) to the four twist field correlation function can now be 
extracted from the OPE 

(T(z)) = + j^d x \ogZ qu (x) + ... (33) 

as z — > x. Evaluating the left hand side by taking the limit z — > x in the expression 

(T(z)) = lim ( g(z,w) - -. 1 ^ = (34) 

x y " ; {z-wfj (z - Xl )(z - x 2 )(z - x 3 )(z - x 4 ) 

1 . / 1 1 1 1 



+ -z/(l - v) + + + 



(z - x ± ) (z - x 2 ) (z - x 3 ) (z - x 4 ) 



we obtain 



Z qu {x) = lim |x 4 | ! " (1 u) (<j v (ti)<j- v (x)<j v (l)<j- v {x4)) = (35) 

const. 

[x{l-x)}^ l -»)[F{x)F{l-x)} 1 / 2 ' 
B. Evaluation of Classical Contribution 

The path integral over the target space fields X % includes a sum over topologically in- 
equivalent configurations from strings wrapping around the compact directions of the torus. 
The main contribution comes from configurations X l cl satisfying the classical equations of 
motion while the effect of fluctuations about these classical configurations is encoded in X qu 
and was calculated in the previous section using conformal field theory techniques. 

In this section we first determine the classical configurations satisfying the equation of 
motion subject to the boundary conditions (4), dictated by the D-brane setup. This is a 
straightforward boundary value problem for the Laplace operator in two dimensions and 
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the solutions can be expressed in terms of holomorphic or antiholomorphic maps from the 
disk onto the target space manifold. We then evaluate the on-shell action and sum over the 
toroidal lattice to obtain the world-sheet instanton contribution to the four point function. 
The solutions to the above boundary value problem are 

dX{z) = awx-^z) = de-^Ui-^z) (36) 
dX{z) = -aui-viz) = -de^u^z) 
dX{z) = bu v (z) = be^-^u^z) 
BX(z) = -bu u {z) = -be^-^Quiz) 

where the coefficients a and b are the only free parameters to be determined. These pa- 
rameters reflect the freedom in specifying the length of the two independent sides of the 
parallelogram. The classical contribution to the path integral is then 

Z c i = e" 5 - (37) 

where 

S d = — *— / d 2 z{dXdX + dXdX) (38) 
Ana' J c+ 

1 7 sin(7rz/)F(a;)F(l - x)(d 2 + b 2 ) 



2na 

where we have used 



Sc. 



d 2 z\Cj v (z)\ 2 = i d 2 z\u 1 _ u {z)\ 2 = 2sin{nu)F{x)F{l -x). (39) 



To determine the coefficients a and b we impose the monodromy conditions 3 

I ds = 2ril ^ Rl [ ds = 2n 2 nR 2 (40) 
J Cl sm(vrz/) Jc 2 

Since X 2 = cot(7rz/)X 1 along C\ 

ds 2 = (dX 1 ) 2 + (dX 2 ) 2 = ( -^-] 2 • (41) 

Similarly, ds 2 = (dX 2 ) 2 along C 2 . A similar calculation as for the quantum monodromy 
conditions then gives 

Cl = 171 [ ; r -— H ; r 42) 

Vsin (nv )F x ) F 1 -x ) V ; 



b = in 



niRi n 2 R 



2 



sin(7rz/)F(a;) F(l — x) 



We assume for simplicity that the branes wrap cycles along the two-torus lattice vectors. 
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Hence, 



2tt 

S d = — sm{nu)F(x)F{l - x) 



riiRi 



sm(iris)F(x) 



+ 



n 2 R 2 



F(l-x) 



The full amplitude is now of the form: 



Z(x)= lim \x i \^ 1 -»\aM°±-Ax)°^)°i-v{x±)) = Z qu V e~ s ^ m ^ 



where Z qu is determined in (36) (up to an overall const, and S c i is defined in (43)). 
Note, in the limit Ri, R 2 — > 00 

e -5 c ;(mi,m 2 ) _^ j 

mi,m2 

and hence the four twist amplitude receives no instanton corrections as expected. 



(43) 



(44) 



(45) 



C. Canonical Form of the Amplitude and Generalizations 

The above calculation was performed for the amplitude of two intersecting branes wrap- 
ping two canonical homology cycles [a] and [b], respectively of the T 2 with the complex 
structure specified by v. We can however reexpress this amplitude in terms of a four-point 
twist amplitude for two branes wrapping two general cycles specified by the wrapping num- 
bers ( rii, mi) and (n 2 ,m 2 ) on T 2 with the trivial complex structure, first, as: 



Z(x) = const.[x(l - x)]- u ^- u \F(x)F(l - x)}' 1 ' 2 



(46) 



Yl ex P sin(7rz/)F(x)F(l - x) 



ri,T2 



\F(x) 



+ 



r 2 L 2 



F(l-x) 



Where Lj are the lengths of the corresponding one-cycles and can be expressed in terms of 
the wrapping numbers and the radii of the torus as: L; L = (niRi) 2 + (miR 2 ) 2 . On the 
other hand sin(7rz/) can be reexpessed in terms of invariant quantities such as the intersection 
number 7 12 = n 1 m 2 — n 2 m x and the lengths L ± , L 2 of the one-cycles as 



sin(7rz/) = 



(2n) 2 I 12 R 1 R 2 
L\L 2 



I12X 



V n i + x 2m i V n 2 + x 2fn 2 



(47) 



where x = R2/R1 is the complex structure modulus. As expected the angle is insensitive 
to the overall scale, and depends only on the wrapping numbers and the complex structure 
modulus. 

It is straightforward to generalize this result to a T 2 with a non-trivial complex structure 
t. In this case it is efficient to parameterize the wrapping numbers in terms of m 8 — > = 
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mi + rrii (see for example [13]). The complete result takes the form (46), but with m/s 
replaced with rfij. 

Of course a generalization of the amplitude to the case of T 2n = T 2 x T 2 ■ ■ ■ (we assume 
the Kahler structure to be a product of the Kahler structures of each T 2 ) is straightforward. 
In this case each twist field is just a product of individual twist fields for each T 2 , and the 
four- twist amplitude is a product of individual twist amplitudes (46). The most interesting 
examples where the above calculations can be applied are cases of Type IIA string theory 
on T 6 = T 2 x T 2 x T 2 with intersecting D6-branes wrapping a product of three one-cycles 
associated with each T 2 . 

For the purpose of performing complete string amplitude calculations it is instructive to 
write down the complete vertex operators for physical bosonic states x an d X* which in the 
(-1) conformal ghost (0) picture 4 : 

3 3 

V-i;x = e " II vi-„e i{1 - Uj)Hj e ik » X " , Vl 1;x * = e"* J] aie^ 1 '^ e ik » x " , (48) 

3=1 3=1 

where Hi corresponds to the bosonized world-sheet fermion ip l (worldsheet super-partner 
of the i-th toroidal coordinate X 1 ). Here we chose to write explicitly the complete vertex 
operator for the bosonic states in four-dimensions; they appear in the internal space at the 
intersection of D6-branes wrapping a product of three one-cycles on T 6 = T 2 x T 2 x T 2 . In 
the case of supersymmetry the intersection angles satisfy the condition [13] Y^=i 7ru i = ^ n 
which ensures that these vertex operators correspond to massless bosonic states, which now 
become superpartners of massless fermionic states with the following (-1/2) super conformal 
ghost picture vertex operators [2]: 

3 3 

V i. =e~ls a f]a{_ u e4-»3) H 3e^ x \ V i. . = A S a ff oje"^ " "Me***" . 

2' ^ , 2' ^ , 

3=1 3=1 

(49) 

Here S a = e 2 1 2 2 and S a = e 2 1 "+" 2 2 represent the spin fields with respective 
positive and negative chirality. [~ e Wl - 2 are bosonized worldsheet fermions ip a with a the 
four-dimensional (complexified) indices.] Note that in the case of supersymmetry the vertex 
operators (48) for x and X* nave the N = 2 worldsheet charge H = X!i=i-^»> +1 an( i — 1> 
respectively and thus correctly represent the vertex operators for the bosonic component of 
the chiral superfield and its complex conjugate, respectively [30]. Similarly, the worldsheet 
charge H for the fermionic vertex operators (49) are respectively —\ and |, again in accor- 



4 Here for clarity we have denoted the bosonic antitwist operators by a\- v modifying the notation a- v 
which was used before for the same operators. 
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dance with N = 2 worldsheet supersymmetry representing the fermionic components of the 
chiral superfield and its complex conjugate, respectively. 

In the above expressions we have suppressed the Chan-Paton factors, however they are 
straightforward to incorporate. The states transform as iV ® M under the U(N) x U (M) 
gauge symmetry of the two intersecting branes (see, e.g., [13] for details). 

The orientifold projection of Type IIA theory involves along with the world-sheet parity 
projection also the mirror symmetry projection, say along the horizontal n-plane of T 2n . 
Note that this projection restricts the value of r to be either or \. Since each brane 
now also has an orientifold image, obtained by a map (rij,mj) — > (rii,—rhi), one can now 
consider four-amplitudes of states appearing at the intersection of a brane, denoted by i, 
and another one, denoted by j*, that is an orientifold image of a brane denoted by j. Note 
that the calculation of the four-amplitudes proceeds analogously as above. It is possible to 
calculate the four-point amplitude of the states appearing at the intersection of, say brane 
i, with its own orientifold image i*. Such states can appear as symmetric or anti-symmetric 
representations of the U(Ni) (for details see [13]). Again, calculation proceeds along the 
same lines. 

In the following we shall determine the crucial normalization const, of the quantum part 
of the correlation function. 

D. Normalization of the Amplitude 

The overall normalization of the four-point amplitude can be determined by factorizing 
the amplitude in the limits x — > or x — > 1 in which the four-point amplitude reduces 
to a product of the two three-point amplitudes. Namely in these limits, the four-point 
amplitude contains a dominant contribution from the exchanges of the intermediate open 
string winding states around the compact directions. In the effective field theory description 
the zero winding states correspond to the exchange of gauge bosons living on the brane 
along the cycle which is not collapsed by the limiting process. As x — > brane 2 contributes 
while as x — > 1 brane 1 contributes. 

For the sake of concreteness we shall focus on a four-dimensional example with D6-branes 
wrapping a product of three-cycles. The physical states at the intersections are represented 
by the vertex operators (48) and (49). For that purpose we shall evaluate the four-point 
disk amplitude S^ki, k 2 , k 3 , fc 4 ) of two bosonic states x an d two x* at the intersections. We 
shall relate this amplitude to the product of two three point functions S 3 (ki, k 2 \ k 3 ) of x an d 
X* states, and the gauge boson via the unitarity condition: 

Qtu u u u\ -r o i2 f d 4k S 3 (k 1 ,k 2 ;k)S 3 (k 3 ,k 4 ;-k) 
S4{ki,k 2 ,k 3 ,k 4 ) =i[g YM \ j ^ _ k 2 + ie ■ ( 50 ) 
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With these preliminaries we now proceed with the calculations of the physical string 
amplitudes. The three-point amplitude takes the form: 

S 3 (h, k 2 ; fc s ) = {V- liX V.^V 0fAlt ) = t^i(2n)^\J2 ^M** " **) " e » ( 51 ) 

V2cr ' — ' 



i=i 

This is the standard three point amplitude, since we have taken into account that 
(cr I/ (0)cri_ l/ (l)) = 1, and that in the picture changing procedure of the gauge-boson ver- 
tex the internal part of the fermionic stress energy tensor does not contribute. We have 
introduced the disk coupling Cm and the coupling g of each vertex operator. The addi- 
tional factor of (2a')~ 1 ^ 2 is due to the picture changing procedure of the gauge field vertex. 
The gauge-field polarization vector is denoted by e M . (For simplicity we calculated the ampli- 
tude only for U(l) gauge field; generalization to U(N) is straightforward.) The factorization 
of the four-point gauge boson amplitude onto the product of two three-point gauge boson 
amplitude yields the standard relationship between Cm and g : Cm = l/(9o a ')- Note that 
this relationship also automatically ensures that on both sides of (50) the dependence on g 
drops out; namely [C D2 gl\ 2 = C D2 gl/a'. 

We have chosen g = e^ 2 which allows one to write the effective kinetic energy action 
for the gauge fields with the pre- factor l/[g Y M} 2 and the kinetic energy for the x fields to 
be canonical. Here g\ M is defined in terms of the full gauge coupling for the specific brane 
as 5 : 

Wym\ 2 = e-*9 2 YM = 2vr J] Sttv^L" 1 (52) 

i=i 

Thus the unitarity condition (50) appears with an extra factor [gy^f] 2 on the right-hand side 
(RHS) of the equation. 

When evaluating the four-point amplitude we chose to picture change the vertex operators 
for both x fields, which in turn ensures that there is no contribution from the internal part of 
the fermionic stress energy contribution. The upshot is the following form of the amplitude: 



S A (h,k 2 ,k 3 ,k A ) = l ^^{2^\Y^k l )Aa' 2 k l ■ h (53) 

i=l 

3 N 

^(l - x)- a ' 1 - 1 Y[[x(l - x)] v ^ 1 - v ' ) Z j (x) + s< — >t 




where the Zj is the four-twist amplitude defined in (46) with v = Vj while s, t are the 
Mandelstam variables. 



See, for example, [38] Vol. II, eq. (13.3.25). 
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In order to compare the LHS and RHS of (50) and thus determine const, we shall evaluate 
the amplitude (46) in the limit x — > first. As x — > 0, F(x) ~ B(u, 1 — v) and F(l — x) ~ 
— log(x/5), where log<5 = 2*0(1) — — — v) and t/j(z) = <ilogr(,2)/<i2. Therefore, to 
take the limit i^Owe must do a Poisson resummation over r 2 . This gives 



Z(x) = ™™J\ W - x)Y^F{x)- 1 (54) 



L2 a/ sin(7rz/) 



exp — 7T- 



mi,m2 



F(l-x) 



— ; sm(7rz/) + 



7rm 2 v / 2o; / 



sin(7rz/) 



2n\fa' const, a/ sin(7rz/) 

^2 A/27T 



E(?) 





r+ 


[ 77fri2 v2a 


'-r\ 











The pre-factors 27rv / a 7 /L 2 from each of the the Zj contribution combine precisely into 
[gy M ] 2 /27r for brane 2 (see eq.(52)). Therefore the contribution of g\ M on both sides of eq. 
(50) cancels. Evaluating the four-point amplitude near x = yields a pole associated with 
s-channel exchange of the corresponding gauge field. Equating the LHS and RHS of (50) in 
turn then determines: 

V2n 



const. 



2*11 



(55) 



3=1 V sin (™j) 

The limit x — > 1 gives a contribution from the t channel exchange of gauge bosons 
associated with brane 1. In this case the resummation is over r\ in (46) which again factorizes 
to [gyjv/] 2 associated with brane 1 in the four- point amplitude (53) and thus cancels the 
same gauge coupling contribution on the RHS of (50). Of course the rest of the calculation 
is consistent with the values of const, in (55). 



E. Generalization of the Lattice Summation 



In the amplitude (46) we assumed that the four twist fields were coming from the same 
intersection and therefore the summation over all possible parallelograms reduced simply 
to a sum over multiples of the lengths, L\ and L 2 , of the two cycles the branes wrap. We 
would like to generalize this amplitude to four twist fields coming from more than just one 
intersection. 

First let us consider the correlation function of a twist-antitwist pair from intersection i 
and a twist-antitwist pair from intersection j. Obviously, the fields coming from the same 
intersection must be separated by a lattice translation and therefore the distance between 
them is again a multiple of the length of one of the two cycles, namely L x or L 2 . However, 
the minimum distance between fields from different intersections is not zero. In particular, 
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it depends on the total number I 12 of intersections between the two branes and the lengths 
of the one-cycles they wrap as we show next. 

By translating the one-cycles by all possible lattice vectors in the covering space C of T 2 
one observes that along one complete cycle each fixed point appears only once and that the 
one-cycle is divided into 1\ 2 equal intervals of length L/I i2 , where L stands for either L\ or 
L 2 depending on the cycle under consideration. Therefore, the minimum distance between 
two different intersections is generally an integer multiple of Lx/I 12 or L 2 /Ii 2 . One must 
first decide on the labelling of the 1\ 2 intersections. There are obviously two evident but 
equivalent labellings, namely, we can index the intersections in increasing order, starting 
from 0, along cycle 1 or cycle 2. Let us, for concreteness, label them along cycle 1. The 
minimum distance between intersections i and j on cycle 1 is then obvious, but the minimum 
distance between the same intersections on cycle 2 is not because the intersection points are 
ordered differently along this cycle. So one must first determine how the fixed points are 
ordered along the second cycle. 

To this end consider the cyclic group of order J 12 , namely 

Z 7l2 = {e, c, c 2 , . . . , c 712 " 1 \c h2 = e} (56) 

where c is the generator of the group. To each fixed point we can uniquely associate an 
element of this group by the rule 

j <— C (57) 

It can then be shown that the ordering of the fixed points along cycle 2 is given by the 
automorphism 

g __> g \ y g e z h2 (58) 

where k is an integer between 1 and I 12 — 1 which depends on the wrapping numbers of 
the two cycles. For this map to be an automorphism obviously k must not divide I 12 for 
then the map is not injective. A general expression for k as a function of the four wrapping 
numbers has proved difficult to find so far apart for special cases of wrapping numbers. 
We emphasize that whatever this expression might be it must guarantee that k does not 
divide I 12 = riim 2 — n 2 m\. Until such an expression is known one can always determine this 
number k by drawing the cycles in the fundamental domain of the torus, k is given by the 
fixed point closest to along cycle 2 (see Fig.4). 

The four-point function with twist fields from intersections i and j gets contributions 
from the two lattice configurations in Figure 5. If di(i,j) oc Lx/I 12 and d 2 (i,j) oc L 2 /I 12 
are the minimal distances between fixed points i and j along cycles 1 and 2 respectively, the 
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FIG. 4: The fundamental domain for a brane with wrapping numbers (3,1) (solid line) and a brane 
with wrapping numbers (1,2) (broken line). There are five intersection points labelled in increasing 
order starting from along the solid brane. Starting from and moving along the second brane 
(broken line) one meets first fixed point 2. This is the integer k that generates the automorphism 
(58) in this example. 



four-point amplitude takes the form 



S A (k u k 2 , k 3 , k A ) = t ^l(2n) 4 5^(J2 • h 



Z) l) (x) + 



( f dxx-*' 8 - 1 ^ - x)-^- 1 T[[x(l - x)] 

v° 

dxx- a,t -\l - x)-*' 3 - 1 Y[[x(l - x)Y^-^zf\x)^ , 



where 



Z (1 \x) = const.[x(l - x)]- u ^ u) [F(x)F(l - x)]- 1/2 



exp sin(TTu)F(x)F(l - x) 



ri,r 2 



riLtY ( d 2 (i,j) +r 2 L 2 



F(x) 



F(l-x) 



and 



(59) 



(60) 



Z^\x) = const.[x(l - x)]- v{l - v) [F(x)F(l - x)]- 1/2 



(61) 
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J 

(a) 



1 

(b) 



FIG. 5: The two configurations for a twist-antitwist pair at intersection i and a twist-antitwist pair 
at intersection j. Both configurations must be included in the string amplitude. 

exp — sin(7ri/)F(x)F(l — x) 

The summation of these two lattice contributions in the case of i=j gives the two terms 
in the amplitude (59). For distinct i and j, however, there is no t-channel massless exchange 
since twist fields from different intersections do not couple. This can be seen explicitly from 
the string amplitude. As x — > or x — > 1 only one of the terms gives a massless exchange 
after Poisson resummation. In particular, for the x — > limit one needs to do a Poisson 
resummation in r 2 to see that only the term survives in this limit. The second term, 
which contains di(i,j), goes to zero even for r\ = in this limit. Analogously, only the 
term contributes in the x — > 1 limit. 

To determine the overall normalization of the amplitude we proceed as in the previous 
subsection. Instead of two s-channel poles and two t-channel poles from gauge bosons living 
respectively on brane 1 and 2, we now get in the amplitude (53) just two s-channel poles, 
one for each type of gauge bosons. The normalization constant is still given by equation 
(55), however. 

One can ask if similar results hold for four-point amplitudes of twist fields coming from 
more than two intersections. Clearly an amplitude of a twist-antitwist pair from intersection 
i with a twist from intersection j and an antitwist from intersection k is not possible since the 
fields coming from intersection i must be separated by a lattice translation which forces k=j. 
However, four-point amplitudes of fields coming from four different intersections are possible. 
In this case there will be a minimum non-zero distance between each pair of twist fields which 
will depend on the particular brane configuration. At most one lattice configuration exists 



/ di(i,j) +r 1 L 1 \ 2 f r 2 L 2 \ 2 
V F(x) J + \F(l-x)J 
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for a given set of twist fields all coming from different intersections. The necessary and 
sufficient condition for a non- vanishing amplitude of two twist fields from intersections i and 
j and two antitwist fields from intersections k and 1 is i-k=l-j. Such amplitudes do not contain 
massless exchanges though and so their overall normalization cannot be determined directly 
by the above method. Nevertheless, this normalization constant is part of the quantum 
amplitude, which is independent of the global effects of the lattice, and hence it must be 
also given by (55). 

III. FOUR AND THREE-POINT FUNCTIONS WITH TWO INDEPENDENT 
ANGLES 

The above method can be directly applied to the problem of a four-point amplitude with 
two independent angles (Fig. 2). The boundary conditions now read 

OX + dX = 0, OX + dX = 0, on (-oo, xi) U (x 2 , x 3 ) U (x 4 , +oo) (62) 
e mv dX + e- inu dX = 0, e~ i7n/ dX + e i7Tl/ dX = 0, on (x u x 2 ) 
e- iwX dX + e inX dX = 0, e inX dX + e~ i7rA <9X = 0, on (x 3 , x 4 ) 

In the appendix we evaluate the quantum amplitude {a^x-^a - v {x 2 )(T -\{x 3 )(T\{x 4 )) . The 
result is 

Z qu (x) = const.x-^ 1 -^! - x)- vX I(x)- 1 ' 2 (63) 

where 

I(x) = (1 - xf- v - x [B(u, X)F 1 (1 - x)K 2 (x) + B(l - u, 1 - A)F 2 (1 - x)^)} .. (64) 

B(v>, A) is the Euler Beta function and F iy Ki are Hypergeometric functions defined in the 
appendix. 

From the boundary conditions (62) we determine the classical solutions 

dX(z) = aw^xiz) ee ae'^-V-,^*) (65) 
dX(z) = -aw^xiz) ee -ae-^-V&^xiz) 
dX(z) = bu u ^ x (z) ee be-™ X Co u ^ x {z) 
§X(z) = -buj u ^ x (z) ee -be™ x u v ^ x {z) 

Again the parameters a and b are arbitrary and reflect the freedom in specifying the 
lengths d\ and d 2 of the four-sided polygon. However, to obtain a three-point amplitude 
one must take the limit x 2 — > x 3 . Unless f — v — A = 0, which is precisely the case of one 
independent angle considered above, one of the two linearly independent solutions becomes 
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singular in this limit. For 1 — v — X > 0, as we will assume without loss of generality, 
u U)1 ^\(z) = [z — xi)~ v (z — x 2 ) u ~ 1 (z — x^) x ~ 1 {z — rc 4 )~ A develops a non-integrable singularity 
at z = X3 in the limit x 2 — ► £3- Therefore, the four-point amplitude that reduces to the 
three point amplitude must have b = 0. This is to be expected since the distance d 2 cannot 
be an independent parameter if one wants to get a three-point amplitude. In fact if b is set 
to zero c?2 becomes a function of d\ and x 2 which tends to zero as x 2 — > £3 as required. To 
keep the discussion general though we first consider the problem with arbitrary a and b. 
The classical action is given by 

1 



Sd — 



Ana' 



~a 2 [ d 2 z\oj 1 ^(z)\ 2 + b 2 I d 2 z\u uA _ x (z)\' 
Kx) 



(66) 



a 2 



;i - x )- 2 ^-% 2 ] 



4a' 

where the integrals have been evaluated using the method of [35] to factorize closed string 
amplitudes into product of open string amplitudes. To determine the coefficients a and b 
we impose the global monodromy conditions 

/ dX 2 = d u I dX 2 = d 2 (67) 

where the contours C[ and C 2 are shown in Figure 3b. The two other lengths of the polygon 
are automatically determined and no other contours provide any additional information 6 . 
Solving these conditions for a and b we obtain 

a = i [(1 - x)- {1 -"- x) B(v, A)Fi(l - x)d x +B(l-u,l- A)F 2 (1 - x)d 2 ] / J{x) (68) 

and 

b = i [B(u, A)Fi(l - x)d 2 + (1 - x) {1 - v ~ x) B(l -u,l- A)F 2 (1 - x)d l ] /J(x) (69) 

where 

= (l- x )-V-»-X) [B(u, X)F l (l - x )] 2 + (l-xf- v -^ [B(l A)F 2 (1 - x)f . (70) 

Although the general four-point amplitude with arbitrary d\ and d 2 is interesting itself, 
we want to evaluate the three-point amplitude which gives the instanton corrections to the 
Yukawa couplings, b then must be set to zero and the monodromy conditions give instead 



6 Note, however, that these two conditions become linearly dependent when 1 — v — A = because d\ = d>2. 
In that case, which is the problem of one independent angle discussed earlier, one should take the contours 
of Fig.3a. 
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while d 2 is x-dependent 



, M B(l-\,l-v)F 2 (l-x) 



= -»>--' B(X, v)F\i\ — x) -° < 72 » 

as a; — > 1 which therefore correctly produces a three-point amplitude in this limit. The 
on-shell action becomes 



4a' ^S(i/,A)Fi(l-a:) 
and the complete four-point amplitude takes the form 

Z 4 (x) = (x>nst.x- v ( 1 - v \l-x)- vX I(x)- 1 ' 2 ^exp-^ ( 

m ^ 



d ^ \ (74) 



B(v,X)F 1 (l-x) 



Here di(m) = L + mL is the most general form the distance d\ can take when the polygon is 
embedded in a lattice. L and the cycle length L will generically depend on the lattice and the 
D-brane configuration. The overall normalization constant can be determined by the same 
method used to fix the normalization of the four-point amplitude with one independent angle. 
We specialize to T 6 since the amplitude of most interest is the one involving factorizable 
three-cycles on T 6 . We find 

const. = 16tt 5/2 . (75) 

The limit x — > 1 of the four-point amplitude produces the full expression for the three- 
point amplitude. In particular, since the conformal weights of the twist fields satisfy h v + 
h\ — h v+ \ = uX, the latter contains the correct singularity in this limit in agreement with 
the operator product expansion 

OiizJOjte) ~ Y. C iP^ z ^ - *i) fcfc ~ fci-fci - (76) 

k 

One can then show that the full three-point amplitude for branes wrapping factorizable 
three-cycles on T 6 takes the form 7 



2 * n 3^)3^,1 < 77) 



J'=l 

3 



J, 



2-na' 



n 167r 2 r(l - Uj)Y{\ - Aj)r(i/j + Xjj v 

j = \ [ r(^.)r(A J )r(i - Vi - a,) ^ 



AAm) 



where Aj(m) is the area of the triangle formed by the three intersecting branes on the j-th 
torus. Note that the amplitude is completely symmetric in all three angles of the triangle 
as it should. 



See Note Added at the end of the paper. 
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The above three-point correlation function of bosonic twist fields (77) is the key contri- 
bution to the physical Yukawa coupling of two fermionic and one bosonic field. The full 
three-point amplitude of these fields can be determined by employing the normalization 
factor for the disk amplitude, Cm = ^/ido®'), and the bosonic string vertex operator in -1 

picture [eq.(48)], go, as determined in Section 2.4. In addition, the normalization factor of 

i 

the fermionic vertex operators in -1/2 picture [eq.(49)] turns out to be 24 yfa' g . [This nor- 
malization can be determined from the string amplitudes of fermionic fields to gauge vector 
bosons, along the same lines as described for the corresponding bosonic fields in Section 2.4.] 
Thus the final expression for the physical Yukawa couplings is given by V2go x Z 3 . 

A comprehensive analysis of the triangles that contribute to this amplitude for a given 
set of intersections in diverse brane configurations and including the effect of non-trivial 
complex structure or Wilson lines has been done in [34]. 

IV. CONCLUSIONS 

We have applied conformal field theory techniques to obtain three-point and four-point 
correlation functions of twist fields from D-branes wrapping factorizable n-cycles of T 2n = 
T 2 x T 2 • • • and intersecting at points in the internal T 2n . The method allows for a complete 
determination of the amplitude including the quantum contribution. Its most interesting 
application is to the three-point function calculation of intersecting D6-branes wrapping 
factorizable three-cycles of T 6 , which in turn gives the complete Yukawa coupling of two 
four-dimensional chiral fermions to the Higgs field. 

The method also applies to the study of the three-point and four-point twist field correla- 
tion functions in models with orientifold and orbifold projections, as discussed in Subsection 
II C. Due to the mirror symmetry projection along an n-plane in T 2n the branes can now 
intersect with an orientifold image of another brane. However the amplitude calculation for 
states at the intersection proceeds analogously. As for the Z 2 x Z 2 orbifold projection, we 
mention that the combined orbifold and orientifold action maps cycles into themselves; thus 
each intersection is accompanied by a combination of the orientifold and orbifold images 
that have to be carefully taken into account. We hope to return to the detailed discussion 
of these contributions in the future work. 

The formalism can be applied to the calculation of three and four fermion couplings of 
classes of Type II orientifold compactifications with intersecting D-branes. Among them 
the four-dimensional Type IIA orientifold compactification with intersecting D6-branes is 
interesting; supersymmetric compactifications of this type are directly related to G 2 com- 
pactification of M-theory. In addition to the Yukawa couplings, the four-point amplitude of 
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chiral fermions can provide low energy corrections to the effective four-fermion coupling and 
may be of phenomenological interest. We plan to address further these effects in concrete 
semi-realistic N=l supersymmetric models such as those constructed in, e.g., [12, 13, 17, 18]. 
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APPENDIX A: EVALUATION OF QUANTUM FOUR-POINT AMPLITUDE 
WITH TWO INDEPENDENT ANGLES 



Here we evaluate the four-point correlation function {a v {xi)a- v {x 2 )a-\{xz)(T\{xi)) . The 
calculation closely parallels the calculation for closed string amplitudes on orbifolds [36]. 
As before we introduce the auxiliary correlators g(z,w), g(z,w), h(z,w) = —g(z,w) and 
h(z,w) = —g(z,w) as defined in section II. A. In terms of 

u v , x (z) = {z- Xl )-»(z - x 2 ) v -\z - x 3 )- x (z - x 4 ) A " 1 (Al) 

these take the form 

( P(z w) \ 
g(z, w) = u^xWw^-xiw) L _ w y + (A2) 

h(z,w) = u v> i-\(z)u Vil -\B(xi) (A3) 

where 

2 

P(z, w)=J2 aijV?* (A4) 

i,j=0 

and the condition 

g(z, w) ~ (z — w)~ 2 (A5) 

determines all coefficients except for a 20 , an and a 02 for which it provides two linear 
equations. Solving these for 002 and an we find 

/rrf \\ v ( I \ 1 ^ A({ Xi }) + z 2 + a 21 z + a 20 

(T(z» = Inn [sM - - ( —^ j = ] —^—^—^—— ) (A6) 

1 s / 1 1 \ 2 1 



+-is(l - v) - - + -A(l - A) . 

v \-v \ / 1-A A 



{z-xi) (z - x 2 ) J \{z - x 3 ) {z-xi) 

where a 2 i = —[(1 — v)x\ + vx 2 + \x 3 + (1 — A)^]. The freedom in a 20 does not affect the 
final result since A({xi}) will be determined after 020 is fixed. Fixing 

a 2 o = ^[(1 - v + X)xiX 3 + (1 - v - X)x x x± - (1 — v - \)x 2 x 3 + (1 + v - X)x 2 xa\ (A7) 
and using SL(2, M) invariance to set x\ = 0, x 2 = x, x 3 — 1, X4 — * 00 we find 

r%ZUx) = i(l-»-A) _ Kl^) _ (1-0(1 -A) (A8) 

where 

= lim -xl x A{{xi}) (A9) 

£4^00 
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and 

Z qu {x) = lim |x 4 | A(1 ^ A) (^(0) ( r_,(x) ( 7„ A (l) ( 7 A ( a ;4)) 

£4^00 

As before A(x) is determined by imposing the quantum monodromy conditions 



/ dX= \ 

Ja Ja 



dX = / dX = 
Jd 

leading again to conditions (22) and (21). We insert 

g(z, w)^[-(l-v- X)x + (l-X)(z-x) + A(x)]a)i_,, A (z) 
g(z,w) -> B(x)uj V;1 - X (z) 

and arrive at the two equations 

[^(1 - v - \)x + A(x)} / dzu)i- Uy \(z) + (1 — A) / dz(z - x)Cui- Vi \{z) 
1 Jd Jd 

-e~ 2 ™B{x) [ dzQx-^z) = 
Jd 

and 

[^(1 - v - \)x + A(x)} / dzu^xiz) + (1 - A) / dz(z - x)lui- v>x (z) 

1 Jd Jd 

-B(x) / dzu x - v ,\{z) = 0. 
Jd 

To solve these we need the integrals 

/ dzu^z) = e^-^Biu, 1 - v)K!(x) 
Jd 

[ dzun-^z) = e^-^Bil - A, 1 - v)(l - x) l ^- x F 2 (x) 
Jd 



where 



K x {x) = F(u,X;l,x) 
K 2 (x) = F(l-i/,l-A;l;x) 
Fi(x) = F(u, A; 1/ + A, x) 



and 



F 2 (x) = F(l - v, 1 - A; 2 - 1/ - A, x). 
One also needs the following identities 

(1 - - X)F(u,X;2,x) = (1 — v — X)K 1 (x) + (1 - x)d x K 1 {x) 
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(1 - - \)F{1 - A, 1 - is; 3 - v - A; 1 - x) = -(2 - v - X)d x F 2 {l - x) (A23) 

B(l-u,l-X)(l-x) 1 - u - x F 2 (l-x) = B(y, X)F 1 (l-x)-B(l-u,l-X)B(u, X)(l-u-X)K 1 (x) 

(A24) 

and 

K^x) = (1 - x) 1 - u - x K 2 {x). (A25) 
After some algebra we arrive at the desired result 

2A(x) = x(l - x)d x log [B(u, A)Fi(l - aO-f^fa) + 5(1 - v, 1 - A)F 2 (1 - x)K 1 (x)} . (A26) 

Hence 

Z qu {x) = const.a^ 1 ^ (1 - x)- uX I( X y 1/2 (A27) 

where 

/(x) = (1 - xf- v - x [B(u, A)Fi(l - x)^^) + 5(1 - i/, 1 - A)F 2 (1 - x)^^)] . (A28) 
For completeness let us also give the result for B(x), namely 



B(x) = \e' 2mX x{l-x) {2 ' v ' x) d x \og 



Bjy, X)F 1 (1 - x)K 2 {x) + B{1 - v, 1 - A)F 2 (1 - x^jx)' 

{K 2 {x)f 

(A29) 
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